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We analyze the conditions under which a statistical error analysis can be carried out in the 
case of TTTT scattering, namely the normality of residuals in the conventional method. Here 

we check that the current and benchmarking analyses only present very small violations of the 
normality requirements. In particular, we show how it is possible to amend slightly the selection 
of the experimental data, and improve the normality of residuals. As an example, we discuss the 
O”*”^ channel and the implications for the /o(500) and /o(980) resonances, obtaining that the new 
selection of data provides very similar and compatible results. In addition, the effect on the /o(500) 
and /o(980) resonance pole parameters is almost negligible, which reinforces the central results and 
the uncertainty analysis performed in these benchmarking determinations. 
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I. INTRODUCTION 

Pion-pion scattering is the simplest hadronic reaction 
in QCD where most theoretical control based on ana- 
lyticity, crossing, chiral symmetry and Regge behavior 
has been achieved (for pedagogical introduction to this 
topic see for example jl|). The most impressive conse¬ 
quence of this development has been the determination 
of S-wave scattering lengths with an order of magnitude 
more precise than the experimental determination m, 
an exceptional case in strong interactions where most of¬ 
ten just the opposite situation is encountered, and exper¬ 
imental precision takes over theoretical predictive power. 

Furthermore, the mass and the width of the lowest res¬ 
onance in the 0++ channel, known as /o(500), and lodged 
in the PDG 0 to stay, have been pinned down unam¬ 
biguously and accurately 0[I1 after a long history full 
of debates. The scalar meson was already proposed by 
Johnson and Teller in 1955 [l^ as the mid-range media¬ 
tor of the nuclear force necessary for nuclear binding of 
atomic nuclei (a historic account on the appearance and 
disappearance of the tr—meson in particle physics may be 
traced from the current 0 and previous editions of the 
PDG booklet). Final values have been reported which 
agree within uncertainties. This high precision has been 
a novel and key convincing culminating element in favor 
of a precise determination of the /o(500) pole parameters. 
Despite of the fact that systematic uncertainties often 
dominate TTTr-scattering analyses, the application of sta¬ 
tistical techniques such as the well-known least squares 
method has been massive and a crucial element to 
quote uncertainties estimates. 
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We remind that the least squares method rests on the 
major assumption that discrepancies between the most 
likely theory and the experiment are statistical fluctua¬ 
tions and more specifically independent random normal 
variables. If there are good reasons to suspect this as¬ 
sumption, one cannot assume the standard approach for 
error propagation. Thus, in order to be able to use 
as a means to obtain uncertainties, this expected nor¬ 
mality of residuals must and can be checked a posteriori. 
Of course, this self-consistency of the fitting procedure 
cannot be decided with absolute certainty, due to the fi¬ 
nite amount of the fitted data. Thus, the key question 
is to decide whether a given finite set of residuals follow, 
within expected statistical fluctuations, a normal distri¬ 
bution pattern or not. The verification of this impor¬ 
tant condition is termed normality test of residuals and, 
though elementary, it has been completely disregarded 
in all previous experimental and theoretical studies in 
TTTT-scattering. This is in contrast with the painstaking 
efforts to reduce errors taken almost elsewhere in these 
works. 

We remind that the theoretical modeling of experimen¬ 
tal data requires within the classical statistical setup an 
assumption both for the signal (the fitting curve) and 
the noise (the statistical fluctuation). However, for a 
self-consistent treatment both assumptions may always 
be checked a posteriori and, if corroborated, error propa¬ 
gation may then be undertaken according to the verified 
distribution of statistical (not necessarily normal) fluc¬ 
tuations. A particular and important example of this is 
provided by the standard and often recommended prac¬ 
tice of re-scaling (by a Birge factor 00) of a too large 
which corresponds to artificially and globally enlarge 
the errors when incompatible data are detected among 
different sets of measurements 00 ). This procedure 
is only justified when we can confidently state that the 
residuals are distributed as scaled gaussians. In such a 
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positive case, we can still propagate errors according to 
the re-scaled gaussian distribution. If, however, this is 
not the case, enlarging the errors in the experiment does 
not entitle to propagate errors according to some distri¬ 
bution. 

In the present paper we want to scrutinize the consis¬ 
tency between tttt scattering data used as the experimen¬ 
tal input information and the theoretical framework used 
to analyze it. This will be done according to the stan¬ 
dard rules of the statistical analysis, namely by inquir¬ 
ing whether or not the difference between the measured 
experimental data and the theory being tested can be 
assumed to be a fluctuation which would likely decrease 
with further and more accurate measurements 0 - 

As it was previously analyzed in Q , there is some ten¬ 
sion among the available experimental tttt data. This 
suggests, as it is usually the case, to make a judicious 
decision on which data should be selected as plausibly 
consistent and, consequently, which data should be dis¬ 
carded as inconsistent. This decision embodies a cer¬ 
tain probability of making an erroneous choice, a num¬ 
ber which can be estimated quantitatively on the basis of 
the assumed model. One naturally expects that the data 
selection improves the statistical behavior of the residu¬ 
als, in which case it seems unlikely that the theoretical 
model contains significant systematic errors. This im¬ 
plies in turn that the uncertainties in the model parame¬ 
ters faithfully reflect the uncertainties of the selected and 
mutually compatible data. 

Our aim here is to restrict the discussion to the simple 
O’*"’'' channel as an example where the necessary statisti¬ 
cal tools can be presented, and to use the Madrid-Krakow 
TTTT analysis [Q-Q as illustration, taking advantage that 
one of us (J.R.E) was one of the authors of this collab¬ 
oration. This procedure can be used as an initial guess 
of a more comprehensive analysis based on solving Roy 
and Roy-like equations supplemented with forward dis¬ 
persion relations and sum rules. As we will show, our 
results based in our restricted analysis will indicate that, 
after the statistical self consistency is over imposed, there 
are no dramatic changes to be expected, providing a sup¬ 
port of the results of Due to the lack of popularity, 

our presentation is intentionally pedagogical 

The paper is organized as follows. In Sect. [II] we review 
the latest tttt analysis performed by the Madrid-Krakow 
group 0. Next, in Sect. imi we review the notion of nor¬ 
mality test as well as some other useful tests. In Sect. IIVI 
we apply these normality tests to the SO wave of the tttt 
analysis in Q. In Sect. El we show how suitably discard¬ 
ing data results regarding normality tests improve. In 
Sect. |Vl] we make a first attempt on making a 3a self- 
consistent fit of mntually compatible data. Finally, in 


^ We note that much used routines such as MINUIT do not 
implement this necessary and simple test and actually very few 
textbooks currently used by particle physicists mention this im¬ 
portant requirement [20|l (one exception is Ref. Ell- 


Sect. Eul we draw our main conclusions and provide an 
outlook for future work. 


II. REVIEW ON THE UNCONSTRAINED FIT 
TO DATA (UFD) 


In a series of works the Madrid-Krakow group 

has used a dispersive approach to built a TTTr-scattering 
amplitude which incorporates analyticity, unitarity and 
crossing symmetry. As a starting point, a set of sim¬ 
ple expressions parametrizing the available experimen¬ 
tal data in each partial wave amplitude were fitted to 
data independently. This is known as Unconstrained 
Fit to Data (UFD). These parametrizations were then 
checked against dispersion relations. When the answer 
was in the affirmative, these parametrizations were used 
as a starting point for a Constrained Fit to Data (CFD), 
in which dispersion relations were imposed as an addi¬ 
tional constraint. Therefore, the set of parameters ob¬ 
tained describe TTTr-scattering consistently with disper¬ 
sion relations. Finally, these parametrizations were then 
used as input for the dispersive Roy and GKPY equa¬ 
tions to perform an analytic extrapolation to the com¬ 
plex plane, yielding precise and model independent de¬ 
terminations of the lightest resonance poles appearing in 
TTTT-scattering [H. 

However, one potential objection to these works, al¬ 
ready noted by the author^ concerns the uncertainty 
treatment. In particular in [Oj, uncertainties were calcu¬ 
lated following two approaches. On the one hand, the 
effect of varying each parameter independently in the 
parametrizations from pi to pi±Api was added in quadra¬ 
ture, disregarding possible correlations among them. The 
errors were further symmetrized to the largest variation 
Api conservatively covering the confidence interval. Al¬ 
ternatively, errors were also estimated using a Monte 
Carlo Gaussian sampling (see e.g. for early im¬ 
plementations of these ideas in tttt scattering) of all UFD 
parameters (within 6 standard deviations). The uncer¬ 
tainties defined by excluding 16% of the upper and lower 
tails of the distribution of a total of 10® events turned out 
be slightly asymmetric. While the essence of the Monte 
Garlo method is to keep correlations, final results were 
not very different from the most naive approach. This 
was partly due to the smallness of errors as well as the 
large number of parameters. 

On the other hand, the way the experimental data are 
represented via standard least squares fitting proce¬ 
dures where one minimizes: 


x"(p) = 


N 

E 

2=1 


Aq^^p 


( 1 ) 


where are experimental observations with estimated 
errors AO™^, Oi{p) is the theoretical model with P pa¬ 
rameters p = {pi,... ,pp) and N the number of data 
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points. At the minimum, 

xLn = minx^(p) =X^(po), (2) 

p 


the most likely theory parameters are po, and the most 
likely theoretical prediction is Of' = Oi(po). The resid¬ 
uals at Po are defined by: 


i?,: = 


OfP - O,(po) 

A^exp 


(3) 


and one req uires the residuals Ri to be normally dis¬ 
tributed (24 [ 2 ^ 0. When normality is fulfilled, the 
test requires (for large N) that within Icr Xmin/^ — 
1 ± \f 2 Jv, with V = N — P the numbers of degrees of free¬ 
dom (dof). Therefore, this requirement precludes both 
loo large and and too small x^/dof. 

In this section we describe the fitting procedure in de¬ 
tail and postpone the normality analysis for the next Sec¬ 
tion. 


A. SO wave parametrization 


Let us first review briefly the different parametriza- 
tions used in Q for the SO wave. The partial-wave is 
written as follows, 


where S^\s) and r]lf'\s) are the phase shift and inelas¬ 
ticity of the partial wave, and k = sjA — is the 
center of mass momentum. 

At low energies, where the elastic approximation is 
valid, a model-independent parametrization ensuring 
elastic unitarity is used. 


AO) _ 
''0 — 




1 


2^ cot(5Q°^(s) — i 


(4) 


In addition, in order to ensure maximal analyticity in the 
complex plane, coti5Q°^(s) is expanded in powers of the 
conformal variable w(s), 


cotSg^^s) = 


2k s— izQ 




Bq + Biw{s) + B 2 w[sY + B 3 w{sy 


j{s) = 


y/s- y/sp - s 
\/s -|- \/So — s 


So = 4M^. 


(5) 


^ These works refer to Nucleon-Nucleon scattering where a set of 
about 8000 neutron-proton and proton-proton scattering data 
below pion production threshold and published from 1950 till 
2013 have been comprehensively analyzed along similar lines 
where many more details can be found. In order to spare the 
effort of the reader unfamiliar with the NN problem we try to 
make the presentation here sufficiently self-contained. 


In the intermediate energy inelastic region, a purely poly¬ 
nomial expansion is used for the phase shift. In addition, 
continuity, and a continuous derivative are imposed at 
the matching point, chosen at s^ = 850 MeV, 


= 


+ \k 2 \{kM — 1 ^ 2 !) > 

(0.85GeV)2 < s <4M|., (6) 

do + B^ + C^+Deis-AMf)-§,, 
4M|- < s < (I.42GeV)2, 


where fc 2 = \/ s/4 — Mj.. Note that we have defined 
kM = |fc 2 (sM), \Sm = S{sm) and 5']^ = dS{sM)/ds, which 
are obtained from Eq. ©■ 

Finally, an elastic SO wave, = I, up to the two- 
kaon threshold is assumed, whereas above that energy, it 
is parametrized as: 


( 0 )/ \ 

Vo io = exp 


-k2{s) 

si/2 


ei + h 


,1/2 


- ^2 

£3 — 

s 


—e 40 (s 


AMf) 


his) 

sl /2 


(7) 


As we will review in the following subsections, the 
number of data points used in Q for each parametriza¬ 
tion is TV = 30,61,28 respectively. In addition, the 
number of parameters of each of them is P = 4, 5,4. 
Thus, it means that to Icr we expect at the minimum 
xVrc = I ± = 1 ± 0-28,1 ± 0.19,1 ± 0.28 with 

V = N — P respectively. Therefore, the total value is 
expected to be jv = 1 ± 0.14. In the following, we will 
review the different sets of data considered in @ to fit 
each of these parametrizations. 


B. Low energy region: conformal fit 

For this fit, it was considered data from the pion 
threshold up to Mk on Kfu decays [1^ , including the lat¬ 
est data from NA48/2 [^, and the selection performed 
in @ of all the existing and conflicting tttt scattering 
data (30l - l3^ between 800 MeV and the matching point 
at 850 MeV (PY05). This selection was obtained from an 
average of the different experimental sets passing a con¬ 
sistency test with Forward Dispersion Relations (FDR), 
where the uncertainties were chosen so that they covered 
the difference between the initial data sets. A detailed 
description of the data selection considered in this aver¬ 
age can be found in @ or in the Appendix of Ref. [s^ ■ 
In Table IXiVl we show the experimental data con¬ 
sidered, together with the fitted values and their cor¬ 
responding residuals. In addition, in Fig. (TJ we show the 
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FIG. 1. UFD for the tttt SO wave. The dark band covers the 
uncertainties, versus the fitted data, namely, the “old” phase 
shift data from decays [2^, the final NA48/2 results [2^ 
and the data average performed in PY05 @]. 

low energy conformal fit in detail, and in Table m we 
show the UFD parameters obtained from the fit. 

C. Intermediate energy region: polynomial fit 

This parametrization describes the SO wave phase shift 
in the energy region between 850 and 1420 MeV. Be¬ 
low the KK threshold, it is fitted to the PY05 average 
of the different experimental solutions compatible with 
FDR [^, and to the re-analysis of the CERN-Munich 
experiments [s^ performed by Kaminski et al. 
which includes a huge estimation of the possible statis¬ 
tical uncertainties. Above the KK threshold, only were 
considered data sets with inelasticities compatible with 
TTTT KK scattering results in the region 4M^ < s < 

(1.25 GeV)^ [S^, namely, the solution (-) of Hyams 

et al. d ata of Grayer et al. [3l| and data of Kamin¬ 
ski et al.js^. Unfortunately, data coming from [s^ only 
provided a statistical uncertainty estimation, so 5° were 
added in 0 as systematic error. In Fig. [21 we show the 
resulting polynomial fit and in Table m the UFD pa¬ 
rameters obtained from the fit. In addition, as we did for 
the previous fit, in Table IXVl we present the experimen¬ 
tal data considered, together with the fitted values and 
their corresponding residuals. 

D. SO inelasticity fit 

For the inelasticity, only data consistent with the dip 
solution were considered, namely the 1973 data of Hyams 
et al. [13, Protopopescu et al. [13 and Kaminski et 
al.[13. However, the data from Kaminski et al. were 
not included in the fit. As explained in [3, the reason 
behind this is that the main source of uncertainty was 
systematic, so the large number of points of Kaminski et 
al. together with the huge statistical errors would lead to 
a fit outcome with much smaller errors than the original 
systematic uncertainties. With only the other two sets. 



FIG. 2. SO UFD, where the dark band covers the uncertain¬ 
ties, versus the fitted data, namely, the data average per¬ 
formed in PY05 [3, the re-analysis of CERN-Munich experi¬ 
ments [13 , and data from [HI . [13 . 


which are incompatible, a fit with a large jd.o.f. was 
obtained, and by rescaling its uncertainties in the inelas¬ 
ticity parameters, it was possible to mimic the dominant 
systematic uncertainties much better. Of course, the re¬ 
sults were still in very good agreement with Kaminski et 
al. In Table llXl of the Appendix, we provide the values 
for the ti parameters, and in Fig. |3|we show the results 
of the unconstrained fit to the SO wave inelasticity data 
up to 1420 MeV. 



FIG. 3. SO inelasticity fit (UFD set) to the tttt —^ tttt scat¬ 
tering data of Hyams et al. (1973) and Protopopescu et al. 
As explained in the text, Kaminski et al. data were not fitted 
[13 although the fit is compatible with them. The dark band 
covers our uncertainties. 
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III. NORMALITY TEST OF THE UFD FIT 
A. General remarks 

In order to understand the discussion below, we have 
to introduce some well known concepts from the statisti¬ 
cal theory of data analysis and specifically on normality 
tests. While the ideas are easy to grasp, these simple but 
extremely informative tests are too often overlooked as 
to deserve to be reviewed in some detail. We avoid on 
purpose the conventional statistical jargon which, after 
exchanging view with other researchers, usually backs 
off theoretical physicists from going in-depth through 
lengthy textbooks and use canned routines as a black 
box. The reader familiar with normality tests may skip 
this section. 

As it is well known, the least squares method implicitly 
assumes that the residuals building the total are nor¬ 
mally distributed. Of course, for a finite number of data, 
we can never be certain about this. We can however, es¬ 
timate the probability by which we would err when we 
decide that this is not the case. In case we have no objec¬ 
tion, we can plausibly assume that further measurements 
will decrease the uncertainties in the fitting parameters; 
otherwise we should expect incompatible determinations 
of fitting parameters. 

In our case we want to decide whether or not the set 
of N residuals (i?i,... Rn) from the x^-fit obey a stan¬ 
dardized normal distribution. The question is whether or 
not the deviations of this empirical sample from the the¬ 
oretical distribution can plausibly be attributed to the 
fact that number of data points N is finite. This can be 
easily done by generating random numbers drawn from a 
normal distribution (.fi,... ,^Ar), € N{0, 1), and com¬ 

paring some quantity of interest T(^i,..., with the 
actually observed empirical result Tobs = T(xi,... ,xn). 
Of course, we can repeat the process as many times 
a = 1,..., M as we want and we will obtain a spread 
of results rth,a = , Cw.a) yielding a distribu¬ 

tion of values for the variable T assuming the variables 
, ^N,a) are uncorrelated normal variables. 

A very important issue is to agree on what quantity 
should be used for this comparison. Let us consider, for 
instance, the moments centered at the origin defined as: 

1 ^ 

flr{xi,...,XN) = (8) 

' 1=1 

Of course, for normally distributed variables 
■ ■ ■ ,^n) becomes a random variable itself which 
follows a certain distribution, 

Pn (/^r) = ) • ■ • ) ^-/v)))) (9) 

where the expectation value is defined as: 

{0)n= / / d^N —... ,^w)- 

J— oo 'y Z'K J— QQ 'y ZtT 

( 10 ) 


This function cannot be computed analytically, except 
for r = 1, 2 or when N ^ oo (featuring the central limit 
theorem 


PNif^r) 


1 r 1^ 

g 2 L J 




( 11 ) 


where the mean (/ir) and the standard deviation 
(A/ir)^ = ifil) — can be expressed in terms of Eu¬ 
ler’s gamma functions [l^. The lowest values r < 8 are 
listed for a quick reference in Table HI Thus, we have that 
if Xi G A[0,1] then within Ict (68%) confidence level and 
for large N, 


1 

— '^Z^ = {Hr)±AHr. ( 12 ) 


When N is not as large (say < 30) one can proceed by 
Monte Carlo by sampling N gaussian numbers and com¬ 
puting the distribution of moments. 


TABLE I. The normalized moments fir of the standardized 
gaussian distribution with their mean standard deviation Afir 
for a sample of size N. For normal distributed data Xi G 
N{0, 1) we expect that to la confidence level = 

fir ± Aflr- 


r 

0 

1 

2 

3 

4 

5 

6 

7 

8 

Mr 

1 

0 

1 

0 

3 

0 

15 

0 

105 

Aflr 

0 

1 

V2 

vTs 

4V6 


3V1130 

3V15015 

240y/^ 

Vn 

y/N 

y/N 

y/N 

y/N 

y/N 

y/N 

y/W 


In the moments method we weight the tails more 
strongly as the power r increases. Therefore, rejecting 
or accepting normality on the basis of each individual 
moment fir tests different features of the distribution. 
This discussion shows that normality is to some extent 
in the eyes of the beholder, and generally it is better to 
use different tests. 

As we have already mentioned, the modeling of data 
requires an assumption both for the signal (the fitting 
curve) and the noise (the statistical fluctuation). In this 
regard, there is an interesting situation where the mo¬ 
ments do not stem from a standardized normal distribu¬ 
tion, but after suitably shifting and re-scaling they do. 
In such a case we can still propagate errors according to 
the shifted and re-scaled distribution. 


B. Normality tests 

The normality test consists of an a priori criterium 
where one decides if the set of data (xi ... ,xn) could 
possibly be normal. To make this decision, one evaluates 
what is the probability p of making an erroneous decision 
of denying the normality of the residuals. The p-value 
is obtained by locating an observable T, known as test 
statistic, from the actual empirical data (a:i ,... ,xi\j) in 
the theoretical distribution that this expected in case the 
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data do follow the normal distribution. A small p-value 
indicates clear deviations from the normal distribution, 
whereas a large p-value indicates that no statistically sig¬ 
nificant discrepancies were found. When comparing Tobs 
to the distribution of T a significance level a is arbitrar¬ 
ily chosen, this determines a critical value Tc- Common 
choices for a are 0.01 and 0.05. We then compare the 
observed T with the critical theoretical value Tc. The 
definition of T in each test determines which inequal¬ 
ity must be fulfilled either from the left or the right. 
Therefore T > Tc 0 the assumption that the empirical 
data xi,... ,a;„ were drawn from the probability distri¬ 
bution p{x) can be rejected with a confidence level of 
100(1 — a)%. If T < Tc there are no statistically signifi¬ 
cant reasons to reject the assumption. This is usually ex¬ 
pressed in short saying that the finite sample {xi,, xn) 
follows the normal distribution iV(0,1) with 100(1 — a)% 
confidence level. 

There are many tests available on the market, and we 
will take here the most popular, which are the Pearson, 
Kolmogorov-Smirnov (KS) and Tail-Sensitive (TS) tests. 
We will also use the moments method already described 
above. 


1. Pearson test 


The Pearson test is one of the first tools used to test 
the goodness of fit by comparing the histograms of the 
empirical data (the residuals of a fit in this case) and of 
the normal distribution. For this purpose a binning must 
be specified and the data is binned accordingly and the 
test statistic is defined as 


Nb 

Tpearson — 

i=l 



/^res _ ^ 

yb 'H ) 

^normal 


(13) 


where N), is the number of bins, nf* is the number of 
residuals on each bin and the expected number 

of data from a normal distribution on the same bin. In 
this case a large value of Tpearson indicates discrepan¬ 
cies between the empirical and normal distribution and 
therefore the null hypothesis is rejected if Tobs > Tc. 

One particular disadvantage of the Pearson test is the 
dependence of Tobs on the number and size of the bins, 
which are set arbitrarily. Different choices can be made 
about the binning, usual methods employ equiprobable 
bins so that is constant; while others use a bin¬ 

ning with a constant bin size (see e.g. Ref. [2^ for more 
details on binning strategies). 


® We are assuming the most frequent case where large T means 
large deviations from the theoretical distribution, as for instance 
in the Pearson and Kolmogorov-Smirnov tests (see below). The 
situation where T < Tc is the relevant inequality can also take 
place, for example in the Tail-sensitive test explained below. 


2. Kolmogorov-Smirnov Test 

The Kolmogorov-Smirnov test (KS) quantifies the dis¬ 
crepancies between two distributions by comparing their 
cumulative distribution functions (CDF). In the case of 
a set of N empirical data Xi, the CDF corresponds to 
the fraction of data that is smaller than a certain value 
X i.e. 

1 ^ 

SNix) =—'^eix - Xi). (14) 

The CDF of the normal distributions is given by 

$(x) = ^ r dte-*^/^. (15) 

V ZtT J — oo 

The test statistic of this test is defined as the maximum 
value of the absolute difference between both cumulative 
distribution functions, that is 

Tks = max |5'Ar(x) - $(x)|. (16) 

— oo<x<oo 

Another advantage of the KS normality test is that a 
good approximation for the p-value exist and is given by 

OO 

fks(Tobs) = 

f=i 

(17) 

This approximation is sufficiently good for > 4 and 
approaches the actual value asymptotically as N becomes 
larger. 

For a normal distribution Xi S A^[0,1] one has that for 
large N and within Icr confidence level, 

Sn{x) = (fix) ± -^\/()){x){l - (f{x)). (18) 

As in the Pearson test, a large value for S'Ar(x) indi¬ 
cates larger discrepancies with the normal distribution 
and therefore the null hypothesis will be rejected if Tks 
is larger than a certain critical values Tc^ks- The critical 
values depend on the sample size N and level of signifi¬ 
cance a with tables for common values of a and < 40 
readily available in the literature. These same tables usu¬ 
ally include a fairly good approximation for A^ > 40, the 
more common case in testing normality of residuals from 
a least squares fit, which we reproduce here as a quick 
reference in table El 


3. QQ plot and tail sensitive test 


The idea of the Tail Sensitive (TS) test comes from the 
normal quantile-quantile (QQ) plot, where one compares 
one to one the empirical points xi < ■ ■ ■ < xn with the 
theoretical points 


n 

(vTi 



dx 






(19) 
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TABLE 11. Large N parametrization of the critical values for 
the Kolmogorov-Smirnov normality test = a/y/N at differ¬ 
ent levels of significance. This approximation is appropriate 
for A < 40 


a 

a 

0.01 

1.63 

0.02 

1.52 

0.05 

1.36 

0.10 

1.22 

0.20 

1.07 


If the empirical points follow a Normal distribution the 
QQ-plot corresponds to the straight line rotated 45°, 
Xn = x^^. A good feature of this representation is that 
the tails are stretched which makes it very suitable to an¬ 
alyze data selection where possible outliers are discarded 
(see Sect. EH)- The TS test defines the test statistic as 

Tts = 2min{min[5i^Ar+i_i(zi), 1 - Bi^N+i-i{zi)\} , 
i 

( 20 ) 

where the cumulative distribution function corresponds 
to the regularized incomplete Beta-function and is de¬ 
fined as 

B..^+i_.(z) = E ( ^ (21) 


and the mapping between 7V[0,1] —>■ 17[0,1] 
1 


Zi = 


J- 


/ 

J —c 


dxe ^ 


( 22 ) 


has been introduced. In Ref. [23 the critical was tab¬ 
ulated for N < 50 and parametrized for iV > 50 as 




(23) 


for the usual significance levels of a = 
0.2,0.1,0.05,0.02,0.01. 


IV. TESTING THE UFD FITS 

After the fitting process and the normality tests have 
been introduced, let us now face the original UFD SO 
wave fit from against them. As mentioned above, this 
is equivalent to check whether the statistical assumptions 
implicit in the standard least squares htting procedure 
are satisfied. Therefore, we will start computing the re¬ 
sulting residuals of the fit, defined by Eq. m where Oi 
denotes the fitted observable, i.e. phase shift or inelas¬ 
ticity data, and that should follow a normal distribution 
within a given confidence level 

As we have pointed out in the previous section, a quan¬ 
titative way to study the normal behavior of the residual 


distribution is to analyze its central moments, which in 
the discrete case are defined through, 

1 ^ 

Mn = ^^(i?*-A„,ea„)”. (24) 

i=l 


In this way, we can compare whether they are compati¬ 
ble within uncertainties with the expected moments of N 
random data points normally distributed, where again, 
N is the total number of residual and the moment er¬ 
ror is given by the square of its variance. We can also 
check whether this situation improves when we consider 
a normalized residual distribution, 


VnR) 


(25) 


where, i?mean = Ri/N is the expected value of the 
residual and V{R) = {Ri — Rrae&nf /{N — I) is its 
variance. In this case, the central moments of the distri¬ 
bution of the normalized residuals are given by: 


I L , 

tin = J^'^{R^- Rmee^nT- (26) 

For completeness, in Table IXVIl we show the experi¬ 
mental data considered, together with the fitted values 
and their corresponding residuals. Since the analysis of 
elastic tttt scattering in the O’*"'*' channel was done sepa¬ 
rately in the low energy region and intermediate region 
and also independently from the inelastic scattering data, 
we have to pose the question on normality in a separate 
fashion. 


A. Low energy region 

In the upper left panel of Fig. HI we show for the confor¬ 
mal low energy fit the distribution of the resulting resid¬ 
uals given in Eq. ([S]) compared with the analytical Gaus¬ 
sian distribution. We can conclude directly from this plot 
that it is not justified to assume that they follow a normal 
distribution. There is an excess of events in the origin. 

In addition, in Table m we show in both cases their 
first moments. By definition, the first two moments are 
/To = I, /J,i = 0. However, for the variance p, 2 , we obtain 
a tiny deviation, which points out that the distribution 
of residuals cannot be considered Gaussian. Since the 
sample of low energy points {N = 30) is rather small, for 
higher values of n, the Gaussian moment uncertainties in¬ 
crease considerably, ensuring that the residual moments 
are going to be satished within uncertainties. 

In the upper right panel of Fig.Hl we show the normal¬ 
ized residual histogram, where we can see, that, despite 
an important improvement, it does not correspond to a 
normal distribution. This can be again checked by com¬ 
paring the last two rows in Table IIIIl where we show the 
central moments of the normalized residual distribution 














FIG. 4. Left: histograms of the resulting residuals after fitting 
the SO wave data. Right: rescaled residual histogram. Upper 
panels: Data described with conformal parametrization. Mid¬ 
dle panels: Data described by the polynomial parametriza¬ 
tion. Lower panels: Inelastic data. The black curve corre¬ 
sponds to an analytical normal distribution. 


and those of a rescaled Gaussian distribution of N ran¬ 
dom points. The first deviations occur for fii. Therefore, 
we can conclude that although only sligh^, the low en¬ 
ergy SO wave phase shift UFD fit from Q violates the 
residual normality test. However, this small violation 
anticipates tiny corrections to the results obtained in Q . 


B. Intermediate energy 

As we did for the low energy case, we will start an¬ 
alyzing the normal behavior of the intermediate energy 
SO wave phase shift fit by studying the behavior of its 
residuals. Their normalized and rescaled distributions 
are plotted in the middle panel of Fig. |4l and point out 
again a excess of residual at the origin, and consequently 
a non-gaussian behavior. As we proceeded for the previ¬ 
ous fit, in order to quantify this deviation, we compute 
again the moments for the normal and rescales residual 
distributions. The results are given in Table HVl 

Since for the intermediate energy region the sample of 
points is bigger {N = 61), the expected errors are smaller 
and the residual constraint becomes this time much more 
stringent. For the non-rescaled case we can see small 
deviation for ^2 and ^ 4 , proving again a non Gaussian 
behavior. For the rescaled case bigger deviations occur 
for n > 3. 



Mo 

Ml 

M2 

M3 

M4 

M5 

Me 

-Rold 

1 

0 

0.6 

0.0 

1.5 

0.1 

4.2 

N(o, Hold 

1 ± 0 

o± 0 

1.0 ± 0.3 

0.0±0.5 2.8±2.1 

0 ± 6 

14 ±30 

-Rnew 

1 

0 

0.9 

-0.1 

2.1 

-0.2 

7 

iV(0.1)„ew 

1 ± 0 

0 ± 0 

1.0 ± 0.3 

0.0± 0.5 2.8 ± 1.9 

-0.1 ± 5 

14 ± 19 


Ao 

Ml 

M2 

M3 

M4 

M5 

Me 

-Rold 

1 

0 

1 

0 

3.5 

0.4 

16 

^(O.l)old 

1 ± 0 

0 ± 0 

1 ± 0 

0.0 ±0.4 2.5 ±0.6 

0 ± 3 

10 ± 7 

^new 

1 

0 

1 

- 0.1 

2.6 

-0.3 

9 

N{0, l)new 

1 ± 0 

0 ± 0 

1 ± 0 

0.0 ±0.4 2.6 ±0.7 

0 ± 3 

11 ± 7 


TABLE III. Central moments of the distribution of residuals 
obtained for the low energy conformal fit versus the moments 
of N random points normally distributed. The subscripts 
“old” and “new” denote respectively the original UFD and 
the UFD for the new selection of data performed in SectionlVl 
fiQ and Hi are by definition 1 and 0 respectively. However for 
fj, 2 , very small deviations appear for the original UFD. For 
higher values of n, they are always compatible due to huge 
uncertainties coming from a small sample of points (N=30). 
In contrast, for the new selection of data all moments are 
compatible with the expected values of a gaussian distribu¬ 
tion. The last four rows contain the central moments of the 
normalized residual distribution obtained in the conformal fit 
versus the moments of N normalized random points normally 
distributed. 



Mo 

Ml 

M2 

M3 

M4 

M5 

Me 

-^old 

1 

0 

0.5 

-0.3 

1.3 

-2 

7 

1^(0, l)old 

1 ± 0 

o± 0 

1.0 ± 0.2 

0.0 ± 0.3 

2.9 ± 1.2 

0 ± 3 

14 ± 12 

.Rnew 

1 

0 

0.8 

0.0 

1.3 

-0.2 

4 

./V(0, Inew) 

1 ± 0 

0 ± 0 

1.0 ± 0.2 

0.0 ± 0.3 

2 . 8 ± 1.2 

-0.2 ± 3 

14 ± 11 


Mo 

Ml 

M2 

M3 

M4 

M5 

Me 

.^old 

1 

0 

1 

-0.8 

5.2 

-12 

53 

^( 0 , lold) 

1 ± 0 

0 ± 0 

1 ± 0 

0.0 ± 0.3 

2.8 ±0.5 

0 ± 3 

13 ± 8 

^^new 

1 

0 

1 

0.1 

2.4 

0.6 

8 

N{0, l)new 

1 ± 0 

0 ± 0 

1 ± 0 

0.0 ± 0.3 

d 

-H 

d 

0 ± 3 

12 ± 8 


TABLE IV. Central non- and normalized moments for the 
residuals of the intermediate energy polynomial parametriza¬ 
tion for the original and new UFD of Section 13 


C. Inelastic data 


In the case of the SO wave inelasticity, the non- and 
normalized distribution of the residuals are plotted in the 
lower panel of Fig.|4l In this case, the main problem does 
not come from an excess of events at the origin, but for 
some values which are far away from it, in contradiction 
with a gaussian behavior. The central moments of both 
distributions are given in Tab. |V] and they show that, 
despite we are working with a small sample {N = 28) and 
the expected errors are big, there are severe deviations, 
some of them of more than an order of magnitude. 
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FIG. 5. (Color online) Rotated Quantile-Quantile plots for residuals of different fits. Left panels: low energy conformal fits. 
Middle panels: intermediate energy polynomial fits. Right panels: high energy inelasticity fits. Upper panels: fit to the original 
database. Lower panels: fit to the newly selected data. The resulting residuals are marked as blue vertical crosses and the 
scaled ones as yellow diagonal crosses. The 95% confidence bands from the TS normality test are also shown (dashed red line). 



Mo 

Ml 

M2 

M3 

M4 

M5 

Me 

-^old 

1 

0 

2.2 

-7.9 

54 

-318 

1981 

A^(0, l)oid 

1 ± 0 

0 ± 0 

1.0 ± 0.2 

0.0 ± 0.4 

3.0± 1.3 

0 ± 3 

14 ± 12 

-Knew 

1 

0 

0.8 

0.0 

1.9 

0.0 

6 

Af(0,l„ew) 

1 ± 0 

0 ± 0 

1.0 ± 0.2 

0.0 ± 0.3 

2 .8± 1.2 

0 ± 3 

14 ± 11 


Mo 

Ml 

M2 

M3 

M4 

M5 

Me 

^old 

1 

0 

1 

-0.8 

5.2 

-12 

53 

N{0, l)oid 

1 ± 0 

0 ± 0 

1 ± 0 

0.0 ± 0.3 

2.8 ±0.5 

0 ± 3 

13 ± 7 

-f^new 

1 

0 

1 

0.1 

2.4 

0.6 

8 

A(0, l)„ew 

1 ± 0 

0 ± 0 

1 ± 0 

0.0 ± 0.3 

2.8 ±0.5 

0 ± 3 

13 ± 7 


TABLE V. Central moments of the residual distribution ob¬ 
tained for the SO wave inelasticity for the original (old) and 
new selection of data (new) (Section 1^. In the case of the 
original UFD, there are very severe deviation with respect to 
a normal distribution. On the contrary, for the new selection 
of data, the agreement is remarkable. 


D. Quantitative checks 

On a more quantitative level, we complement the pre¬ 
vious results with Tables IVII VIIIVlUl for the Pearson, KS 
and TS tests respectively. It is interesting to note that 
for the Pearson and KS tests the p-value is very high for 
the Inelastic fit, whereas the TS gives a too low value. 
This is mainly due to the outlier with Ri = —6.24 sit¬ 
ting at the tail of the distribution. We also note that for 
these small data points values (N ), the power of the test 
(probability of not giving a false positive) is low. Finally, 
QQ-plots are shown in Fig. [S] All this reinforces the lack 
of normality of residuals. 


TABLE VI. Results of the Pearson normality test of the resid¬ 
uals obtained by fitting the different data sets with the confor¬ 
mal, polynomial and inelastic parametrizations. The results 
of the test of the scaled residuals for every case is shown below 
the corresponding line. The critical value Tc corresponds to 
a significance level of a = 0.05. 


Database 

Fit 

N 

Tc 

Fobs 

p-value 

Conformal 

UFD 

30 

14.07 

14.80 

0.039 





10.43 

0.021 

Conformal new 

UFD 

25 

14.07 

5.40 

0.611 





5.40 

0.611 

Polynomial 

UFD 

61 

18.31 

21.05 

0.021 





6.98 

0.727 

Polynomial new 

UFD 

52 

16.92 

4.15 

0.901 





3.00 

0.964 

Inelastic 

UFD 

28 

14.07 

5.14 

0.642 





6.71 

0.573 

Inelastic new 

UFD 

27 

14.07 

2.93 

0.892 





2.93 

0.892 


V. NEW SELECTION OF DATA 

The results presented in Section Hvl highlight that the 
UFD fits to the tttt scattering amplitude in the O'*''*' chan¬ 
nel from Q do not satisfy a posteriori the statistical as¬ 
sumptions implicit in the minimization. This incon¬ 
sistency in the fit can be due to either the presence of 
mutually incompatible data or a biased choice of the fit- 
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TABLE VII. Same as table IVTl for the Kolmogorov-Smirnov 
test 


Database 

Fit 

N 

Te 

Fobs 

p-value 

Conformal 

UFD 

30 

0.242 

0.253 

0.036 


UFD 



0.216 

0.105 

Conformal new 

UFD 

25 

0.264 

0.158 

0.513 


UFD 



0.157 

0.520 

Polynomial 

UFD 

61 

0.174 

0.170 

0.053 


UFD 



0.097 

0.582 

Polynomial new 

UFD 

52 

0.189 

0.082 

0.846 


UFD 



0.058 

0.991 

Inelastic 

UFD 

28 

0.250 

0.118 

0.785 


UFD 



0.186 

0.254 

Inelastic new 

UFD 

27 

0.254 

0.140 

0.611 


UFD 



0.129 

0.715 


TABLE VIII. Same as table IVll for the tail sensitive test 


Database 

Fit 

N 

Tc 

Fobs 

p-value 

Conformal 

UFD 

30 

0.0039 

0.0049 

0.061 


UFD 



0.0155 

0.164 

Conformal new 

UFD 

25 

0.0042 

0.1662 

0.794 


UFD 



0.1686 

0.799 

Polynomial 

UFD 

61 

0.0027 

0.0010 

0.020 


UFD 



0.0081 

0.132 

Polynomial new 

UFD 

52 

0.0029 

0.0912 

0.681 


UFD 



0.3101 

0.988 

Inelastic 

UFD 

28 

0.0040 

1.2 X 10"® 

0.0002 


UFD 



0.001 

0.016 

Inelastic new 

UFD 

27 

0.0040 

0.1505 

0.769 


UFD 



0.1913 

0.851 


ting curve. We will explore the consequences that a new 
selection of data implies regarding normality of residuals. 


A. New selection of low energy data 

As we have seen before, there are three different sets 
of experiments included in the low energy region. The 
NA48/2 set corresponds to the newest data on K 14 , 
decays and incorporates a rigorous statistical and sys¬ 
tematic error analysis, so there is no reason to modify 
or alter this set of data. As we have commented before, 
the PY05 average selection corresponds to an average of 
the different experimental solutions that passed a con¬ 
sistency test with FDR and other sum rules. The large 
uncertainties assigned covered the difference between the 
different data sets, and at some point, they can be con¬ 
sidered arbitrary. Since, as we can see in upper Fig. 21 
there is an excess of residuals at the origin, it is reason¬ 
able to assume that they may be overestimated. The last 



s‘“ (MeV) 

FIG. 6. Old versus new UFD at low energies, the dark band 
covers the nncertainties of the old parametrization. 


set of data considered in the low energy fit corresponds 
to old data on K{,i decays from [ 2 ^ , whose precision can¬ 
not be consider as the same level than the recent NA48/2 
analysis. Therefore, it makes sense to get rid of those old 
data points on decays which lie within uncertainties 
in the same energy bin than a NA48/2 value. Thus, in 
order to improve the Gaussian check of the fit, we divide 
by two the uncertainty of the PY05 data average and 
get rid of the old K decay corresponding to the following 
values of {289, 317, 324, 340, 367}. The comparison 
of the tiny differences between the old and the new fit is 
plotted in Fig. [51 Its corresponding parameters are given 
in Table HXl 

Finally, in the upper panel of Fig. 0 we plot the resid¬ 
ual distribution for the new fit. Despite there are still 
some deviations, there is a clear improvement, which can 
be again checked by computing the central moments of 
the distribution. They are given in Table uni showing 
a complete agreement between them and the expected 
values of N random points normally distributed. 

In this was, by imposing the Gaussian check in the 
data selection, we are able to constraint part of the ar¬ 
bitrariness assumed in the data selection, namely, the 
systematic uncertainty taken for the PY05 average and 
the old K decay data points considered. 


B. New selection of intermediate energy data 

On the one hand, as we can can see in middle Fig. 21 
there is again a excess of residuals in the central region, 
i.e. for i? ~ 0, which could be considered as a signal 
of overestimated uncertainties. On the other hand, the 
uncertainties of the PY05 average data, the systematic 
error added to Hyams and the errors given by Kamin¬ 
ski et al. can be considered at some point arbitrary. 
Therefore, in order to improve the normal behavior of 
the residual distribution, we reduce by a factor 1/2 the 
uncertainties of the PY05 average, we sum 3° instead 
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of 5° as systematic uncertainty to the data from Hyams, 
and rescale by a factor of 3/4 the uncertainties of Kamin¬ 
ski data. In addition, since there are several Kaminski 
data points which are compatible with those of Grayer, 
Hyams or PY05 but with higher uncertainties, we get 
rid of the Kaminski data point at the following ener¬ 
gies {930,970,1010,1050,1210,1390}. The normalized and 
rescaled distribution of the residuals for the new data se¬ 
lection are plotted in middle Fig. [3 whereas the central 
moments are given in Table IIVI Despite there are still 
some deviations, the improvement is clear. The value of 
the parameters corresponding to the new polynomial fit 
are given in Table m and the comparison between the 
new and old polynomial curve is plotted in Fig. [H 


C. New selection of inelasticity data 


In this case, as we can see in the histograms depicted 
in the lower panel of Fig. |4l rather than by overestimated 
errors, the problem comes from incompatible data points, 
and in particular for the Hyams data point at y/s = 994 
MeV. If we simply get rid of this point, the new fit im¬ 
proves striking the normal behavior as we can see in the 
lower panel of Fig. [71 The central moments of the new fit 
are given in Table IVl and show a perfect agreement. The 
new inelasticity parameters are given in Table m and 
the difference between the old and new parametrizations 
are plotted in Fig. jS) 



FIG. 8. Old versus new UFD for the polynomial parametriza- 
tion, the dark band covers the uncertainties of the old 
parametrization. 



FIG. 9. Old versus new UFD for the inelasticity, the dark 
band covers the uncertainties of the old parametrization. 


D. Resonance pole parameters 


As it has been pointed out in the introduction, one 
of the main achievements of the dispersive tttt analy¬ 
ses based on Roy and Roy-like equations H, 0, 0 has 
been the precise and unbiased determination of the low¬ 
est scalar resonance states ELIil, namely, the /o(500) 
and/o(980). In particular, in HI , the CFD parametriza¬ 
tions of 0 were used as input for the analytic extrapola¬ 
tion to the complex plane of the dispersive once- and 
twice-subtracted Roy equation, and thus, allowing to 
look for the lowest-lying poles on the second Riemann 
sheet. Therefore, it is relevant to check the effect of the 
new selection of data on the scalar pole determinations. 
In Table 1 X 1 we show the pole positions of the /o(500) and 
/o(980) resonances for the old and new UFD, obtained 
from the twice-subtracted Roy equations. For complete¬ 
ness, we also show the results for the p{770) resonance. 
Note that the values of the p pole depend on the new SO 
wave data selection through the dispersive integral. The 
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SO wave 

UFD 

New UFD 

Bo 

7.26 ± 0.23 

7.38 ± 0.23 

Si 

-25.3 ± 0.5 

-24.5 ± 0.5 

B2 

-33.1 ± 1.2 

-35.4 ± 1.2 

Bo 

-26.6 ± 2.3 

-34.7 ± 2.7 

zq 


M, 


0.645 

1.111 

1 - 

0.723 

0.691 

1 + vWl' 

1.277 

1.309 

do 

(227.1 ± 1.3)° 

(228.8 ± 1.3)° 

c 

(-660 ± 290)° 

(-466 ± 294)° 

B 

(94.0 ± 2.3)° 

(85.1 ± 2.3)° 

c 

(40.4 ±2.9)° 

(60.7± 2.9)° 

D 

(-86.9 ± 4.0)° 

(-92.3 ±4.0)° 


0.552 

0.761 

1 - -splfd 

0.811 

0.794 

1 -f •jtjT’ 

1.189 

1.206 

€l 

4.7 ±0.2 

4.6 ± 0.1 

€2 

-15.0± 0.8 

-14.6 ± 0.4 

€3 

4.7 ±2.6 

3.9 ±4.4 

e4 

0.38 ± 0.34 

0.37 ±0.20 


2.668 

1.039 

1 - vWl- 

0.711 

0.705 

1 -f sfifd 

1.289 

1.294 


TABLE IX. SO wave parameters for the old and new UFD 
sets, indicating also the number of fitted data N, the cor¬ 
responding expected value boundaries 1 ± yj%jv. 

The first set corresponds to the low energy parametrization, 
\fs < 0.85 GeV, the second set to the intermediate energy 
conformal fit and the last one to the inelastic parametriza¬ 
tion up to \fs = 1.42 GeV. 

small differences between both determinations point out 
the small effects of the new data selection and reinforce 
the results obtained in [H. 


old UFD(MeV) new UFD (MeV) 

^ 5j-p(500) (482 ± 16) - i(268 ± 16) (483 ± 17) - i(266 ± 16) 

^/o(980) (loot ±3)-i(8 ±9) (1001 ± 5) - i(12 ± 10) 

\/p(770) (764.4 ± 2.0)-i(73.3± 1.2) (764.2 ± 1.9) - i(73.2 ± 1.2) 

TABLE X. /o(500), /o(980) and p(770) resonance pole po¬ 
sitions of the old and new UFD obtained from the analytic 
continuation to the complex plane of the twice subtracted 
Roy equation . The tiny differences between both determi¬ 
nations strengthen the error analysis carried out in [^. 


VI. ELEMENTS OF A SELF-CONSISTENT FIT 

As we have mentioned and demonstrated in the previ¬ 
ous section, the selection of data is an essential ingredient 
for a self-consistent fit, namely a fit where residuals fol¬ 
low a normal distribution. However, this selection has 
been done by hand, and a more general procedure would 
be most useful. Guided by previous experience in NN 
scattering we try to explore a similar selection 

method. 


Of course, one important aspect is that, unless specifi¬ 
cally proven incorrect, all available data as they are pub¬ 
lished by the experimentalists should be taken into ac¬ 
count to make the data selection. However let us re¬ 
member that, in the particular case of tttt scattering, 
phase shift and inelasticity data are not obtained from 
direct measurements, but indirect, requiring the use of 
models, and thus, involve systematic uncertainties. Note 
by instance, that the PY05 average 0 used in both, 
the conformal and intermediate-energy parametrizations 
of 0 was obtained confronting these experimental anal¬ 
yses against Forward dispersion relations. 

Nevertheless, if in general we assume from the start 
that all experiments are correct; our goal is to reach a 
consensus among published experimental data which are 
or could be mutually compatible. This is an important 
assumption, but it is the most objective one we can make 
without any detailed knowledge on how the experimen¬ 
tal analysis was carried out; if we did know, we would 
decide based on our own judgement. We hope, however, 
that given a sufficiently large number of data, statistical 
analysis will help us to make a judicious choice and to 
reach a consensus among individual data measurements. 
In this way, we envisage the possibility that simple ex¬ 
perimental values in a given experiment can individually 
be discarded without questioning the whole experiment. 
A practical way to implement this method is to follow 
the 3cr self-consistent procedure proposed in [2^, which 
is given by the following three-steps algorithm: i) Make a 
first and global fit, ii) Discard from all data those fulfill¬ 
ing Rf > 9 iii) Re-fit the remaining data and go to step 
ii) until convergence. This selection method will gener¬ 
ate a boundary between accepted and rejected data and 
there will be a flow of data across the boundary during 
the iterative process. 

The good feature of this method is that in the case of 
two mutually incompatible data it helps to decide which 
one is better suited even when initially and individu¬ 
ally one would reject both. Of course, this decision is 
controlled by the fitting theory and unforeseen restric¬ 
tions on the theory may induce an undesirable bias in 
our choice. It is thus important that the fitting curve 
is flexible enough to prevent this situation @. Note that 
after this process has been carried out, we have still to 
check whether the residuals pass the normality test. 

In order to check the disagreement between the tttt 
data sets considered in the previous sections, we ap¬ 
ply the self-consistent approach to the SO wave phase 
shift, performing a global fit to the data sets described in 
Sect II V Al and IIV Hi The advantage of a global fit, instead 


From this point of view it is preferable to have an excess of fitting 
parameters, since their redundancy will emerge through correla¬ 
tions among them. In the opposite situation, a too restrictive 
choice may reject data just on a lack of flexibility. Note that we 
want to use the most general and admissible theory which can 
accommodate all possible but correct data. 
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of considering separately the low and energy region, is 
that we can treat on an equal footing both energy regions, 
and thus, without introducing a priory any bias into our 
analysis. Note, however, that the hierarchy imposed in 
the original work of Q, allows to take advantage of the 
high precision accomplished in the latest NA48/2 [ 2 ^. 

In all, we find that if we fit the data in the SO up 
to y/s = 1.42 GeV convergence is achieved already with 
two iterations, with a total of only three points b eing 
discarded, namely, the “old” K decay data point 
at y/s = 381.4 MeV, (Table IXIVj) and the Kaminski 
et al. [13 data points at ^/s = {970,1050} MeV, (Ta¬ 
ble |XV]). In Table IXlIIl we provide the resulting param¬ 
eters in the first and second iterations. Note that the 
differences between the new UFD parameters and the 
original ones Table IIXI are due to the global fit we are 
performing in the self-consistent method. 

As compared to the more detailed discussion in the 
previous section, it is fair say that the 3cr self-consistent 
approach is much simpler and effective, and does not re¬ 
quire from the theoretician’s side a detailed discussion 
of the published data but rests on the assumption that 
most of data are correct 0 and that statistical regularity 
does imply a consensus among the mutually compatible 
data. This makes sense of course provided the theoretical 
model used to undertake the analysis is flexible enough. 

Since the selection of data implies a sharp cut for resid¬ 
uals with i?? > 9, a comparison with the standard nor¬ 
mal distribution is not appropriate; the empirical dis¬ 
tribution will have no tails. However it is possible to 
compare against a truncated normal distribution iV(0,1) 
with V(a:) = 0 for |x| > 3. A first comparison between 
this truncated distribution and the resulting residuals has 
been made by looking at the corresponding moments. Of 
course the values in Table |T] are no longer valid due to the 
subtraction of the tails. The moments fir and A fir for 
a finite size sample of the 3cr-truncated distribution are 
shown in Table IXll The results for N = 88 along with the 
empirical values for the residuals are shown in table IXIII 


r 

0 

1 

2 

3 

4 

5 

6 

fl-r 

1. 

0. 

0.973 

0. 

2.680 

0. 

11.240 

A/Xr 

0 

0.986 

1.316 

3.352 

7.215 

18.928 

47.331 

\/N 

v/N 

\/N 

v/N 

v/N 

v/N 


TABLE XI. The normalized moments fir of the standard¬ 
ized 3(j-truncated gaussian distribution with their mean stan¬ 
dard deviation Afir for a sample of size N. For normal dis¬ 
tributed data Xi € iV(0,1) we expect that to la confidence 
level /-V = /ir ± Aflr- 

It is notable the improbably low jv values for this 
new fit shown in table IXIIII This may be a consequence 


® However, we stress again that in the case of 7r7r-scattering, where 
systematics dominate the experiment analyses, dispersion rela¬ 
tions allow to discard data sets which are in contradiction with 
the dispersive constraints Q. 


_ na Ui _ U 2 _Ms_ fM _Ms_ Me 

1 -0.11 0.59 -0.58 1.75 -3.41 9.36 

iV(0,l) 1±0 OibO.ll 0.97±0.14 0.0± 0.36 2.7±0.8 0±2.1 11.24±5 

TABLE XII. Central moments of the distribution of residuals 
obtained for the global self-consistent fit versus the moments 
of A = 88 random points normally distributed with trunca¬ 
tion at |x| = 3. fio and fii are by definition 1 and 0 respec¬ 
tively. However from fi 2 to fi 5 , deviations appear for empirical 
distribution. The deviations are larger for odd moments given 
the large asymmetries of the residuals. For n = 6, the values 
are compatible due to higher uncertainties. 

of only discarding data with a too large contribution to 
the (mostly underestimated errors). Data with over¬ 
estimated error bars remain in the fit and their too low 
contribution to the pulls down the final value. This 
actually shows that a selection of data does not necessar¬ 
ily complies with (truncated) normality. 


SO wave 

Fitl 

Fit2 

Bo 

6.32 ± 0.28 

6.10 ± 0.27 

Bi 

-21.6±0.6 

-20.1 ± 0.6 

S2 

-32.4 ± 1.5 

-31.9 ± 1.4 

Bo 

-43.9 ±3.3 

-50.8 ± 3.2 

-20 



do 

(227.2 ± 1.3)° 

(228.8 ± 1.3)° 

c 

(-675 ± 248)° 

(-195 ± 248)° 

B 

(93.7 ± 2.3)° 

(92.2 ± 2.3)° 

G 

(49.4 ± 2.9)° 

(52.2 ± 2.9)° 

D 

(-86.4 ±4.0)° 

(-89.8 ± 4.0)° 

N 

88 

85 

^ jv 

0.54 

0.38 


TABLE XIII. SO wave parameters for the self consistent global 
fit. Note that in this case, the whole phase shift energy region 
has been fitted at the same time. Fit 1 and Fit 2 stands 
for first and second iterations, indicating also the number of 
fitted data N and the corresponding jv (expected value is 
1±^2 /m. Convergence is achieved in the second iteration.The 
first four lines correspond to the low energy parametrization, 
v/s < 0.85 GeV, and the last nine to the parametrization up 
to yi = 1.42 GeV. 


VII. CONCLUSIONS 

One of the most relevant accomplishments in hadronic 
physics in the last decade has been the reliable deter¬ 
mination of the mass and width of the lowest scalar 0++ 
resonance, also known as the cr-meson, which has entered 
into the PDG as the /o(500) state. This has occurred 
as a side-product of comprehensive long-term studies in 
TTTT scattering, a particularly simple reaction where many 
theoretical constraints such as crossing, analyticity, uni- 
tarity, chiral symmetry and Regge behavior allow for con¬ 
vincing and accurate theoretical predictions. This con¬ 
clusion holds after a long tour the force, and thus the 
significance of this major issue should not be underesti¬ 
mated. In the present work we have re-analyzed the sta- 
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tistical treatment of tttt scattering data from the point 
of view of normality tests; which have so far been over¬ 
looked. The basic aim was to check whether the currently 
accepted and agreed bench-marking analyses carried out 
during the last decade fulfill the condition that the dif¬ 
ference between the fitted data and the fitting theoretical 
curve can be regarded as a statistical fluctuation. This is 
an elementary requirement on the applicability of statis¬ 
tical methods such as the least squares fit which can only 
be carried out a posteriori. As an example, we have ap¬ 
plied several conventional normality tests to the SO wave 
of the precise tttt analysis performed in Q. This test 
has pointed out a tiny violation of the normality require¬ 
ments. When the normality test fails many questions 
should be asked, but the most obvious ones address ei¬ 
ther the compatibility of the data base used to carry out 
the x^-fit or the incapability of the theory to describe the 
data or both. We have carried out a preliminary analysis 
along these lines. While our study can definitely be im¬ 
proved, we have analyzed several strategies incorporating 
many of the elements that a full analysis might contain, 
including data selection and normality tests. However, 


the small differences obtained suggests that there is no 
need for a new reanalysis of tttt scattering. 

We find that by changing slightly the selection of 
data, normality of residuals can be achieved in a sig¬ 
nificant way. This allows to propagate errors and hence 
to reassess the estimation of uncertainties in the scalar 
resonance parameters. We find small and compatible 
changes, reinforcing ultimately the benchmarking deter¬ 
minations carried out during the last decade, and in par¬ 
ticular the results performed in Q. 
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TABLE XIV. Experimental SO phase shift data used for the 
low-energy conformal fit, together with the fitted values and 
the residuals obtained from Eq. [3] The experimental data 
can be gathered in three different blocks. The newest Ke 4 
data from NA48/2 [2^, which incorporates a rigorous sys¬ 
tematic and statistical error analysis, Ke 4 data from old ex¬ 
periments [ 2 ^, and the average result collected in PY05 [^. 


Appendix A: Data and residuals 


Experiment 


(O^P ± AO®’''’)'"' 

^th ( 0 ) 

Ri 

NA48/2 

286.1 

2.32 ± 1.89 

2.91 

-0.31 


295.1 

4.67 ± 1.26 

4.83 

-0.12 


304.9 

6.20 ± 0.95 

6.30 

-0.11 


313.5 

7.57 ±0.82 

7.65 

-0.10 


322.0 

8.70 ± 0.64 

8.96 

-0.41 


330.8 

9.99 ± 0.58 

8.96 

-0.55 


340.2 

12.44 ± 0.75 

11.75 

0.91 


350.9 

13.69 ± 0.57 

13.44 

0.44 


364.6 

16.60 ± 0.55 

15.62 

1.78 


389.9 

19.43 ± 0.50 

19.83 

-0.80 

K-decays 

285.2 

-1.63 ± 2.28 

2.73 

-1.91 


289.0 

3.36 ± 7.44 

3.44 

-0.02 


299.5 

3.36 ± 1.42 

5.43 

-1.46 


303.0 

11.53 ± 4.00 

6.00 

1.38 


311.2 

7.25 ± 1.08 

7.31 

-0.06 


317.0 

7.08 ± 2.85 

8.19 

-0.39 


324.0 

8.87 ±0.96 

9.27 

-0.41 


335.0 

11.29 ± 2.28 

10.95 

0.15 


340.4 

12.05 ± 0.84 

11.80 

0.29 


367.0 

15.80 ± 2.77 

16.02 

-0.09 


381.4 

16.86 ± 1.18 

18.40 

-1.31 


810 

88 ± 6 

86.36 

0.27 

830 

92 ± 7 

89.20 

0.27 

850 

94 ± 6 

92.19 

0.27 

870 

91 ± 9 

95.5 

-0.50 

910 

99 ± 6 

104.0 

-0.83 

912 

103 ± 8 

104.6 

-0.19 

929 

112.5 ± 13.0 

110.1 

0.18 

935 

109 ±8 

112.6 

-0.45 

952 

126 ± 16 

121.8 

0.26 


PY05 






16 


Experiment 


(0“P ± 

^th (o) 

Ri 

PY05 

870 

91 ± 9 

95.5 

-0.50 


910 

99 ±6 

104.0 

-0.83 


912 

103 ± 8 

104.6 

-0.19 


929 

112.5 ± 13.0 

110.1 

0.18 


935 

109 ± 8 

112.6 

-0.45 


952 

126 ± 16 

121.8 

0.26 


965 

134 ± 17 

132.7 

0.10 


970 

141 ± 21 

138.4 

0.14 

Kaminski 

910 

101.0 ± 9.1 

104.0 

-0.22 


930 

113.2 ± 7.2 

110.5 

0.25 


950 

123.2 ± 7.1 

120.5 

0.19 


970 

110.1 ± 10.4 

138.4 

-1.81 


990 

205.2 ± 7.5 

192.6 

1.11 


1010 

233.5 ± 21.3 

230.6 

0.09 


1030 

178.5 ± 47.9 

234.7 

-0.78 


1050 

192.2 ± 11.5 

239.1 

-2.72 


1070 

257.0 ± 11.9 

243.8 

0.74 


1090 

246.7 ± 20.3 

248.7 

-0.06 


1110 

262.4 ±8.9 

251.4 

0.81 


1130 

221.8 ± 22.9 

253.8 

-0.93 


1150 

267.4 ± 16.1 

256.3 

0.46 


1170 

274.6 ± 18.2 

259.1 

0.57 


1190 

258.0 ± 21.9 

262.2 

-0.13 


1210 

265.2 ± 13.5 

265.5 

-0.02 


1230 

255.7 ± 11.5 

269.2 

-0.78 


1250 

258.2 ± 14.8 

273.1 

-0.67 


1270 

258.5 ± 11.7 

277.4 

-1.08 


1290 

273.2 ± 10.9 

282.0 

-0.53 


1310 

297.3 ± 9.3 

286.9 

0.74 


1330 

279.6 ± 16.1 

292.2 

-0.52 


1350 

305.6 ± 10.3 

297.9 

0.51 


1370 

296.9 ± 21.7 

303.9 

-0.32 


1390 

321.8 ± 17.8 

310.4 

0.13 


1410 

321.2 ± 14.2 

317.3 

-0.36 

Hyams 

1020 

230.2 ± 6.4 

232.6 

-0.38 


1060 

242.1 ± 6.4 

241.4 

0.10 


1100 

246.6 ± 6.2 

250.5 

-0.63 


1140 

255.3 ± 6.1 

255.0 

0.04 


1180 

256.3 ± 5.5 

260.6 

-0.78 


1220 

264.9 ± 5.3 

267.3 

-0.45 


1260 

274.1 ± 5.8 

275.2 

-0.19 


1300 

283.9 ± 5.5 

284.4 

-0.09 


1340 

289.1 ± 5.3 

295.0 

-1.11 


1380 

309.6 ± 5.2 

307.1 

0.47 


1420 

322.2 ± 6.7 

320.9 

0.17 

Grayer 

991.7 

231.3 ± 13.5 

212.8 

1.35 


1002.3 

206.6 ± 20.1 

229.1 

-1.11 


1031.7 

225.0 ± 22.5 

235.1 

-0.45 


1051.0 

245.3 ± 12.1 

239.4 

0.49 


1070.8 

249.8 ± 10.3 

244.0 

0.56 


1091.7 

251.6 ± 7.7 

249.1 

0.31 


1112.5 

251.6 ± 10.3 

251.8 

-0.03 


1133.3 

257.9 ± 12.1 

254.2 

0.30 


1150.0 

259.6 ± 12.1 

256.3 

0.27 


1173.3 

261.9 ± 10.3 

259.6 

0.22 


1193.3 

259.6 ± 8.1 

262.7 

-0.38 


1208.0 

270.0 ± 6.0 

265.2 

0.81 


1225.6 

271.0 ± 7.4 

268.3 

0.36 


1246.4 

282.3 ± 6.5 

272.4 

1.53 


1264.4 

278.2 ± 8.6 

276.2 

0.23 


1290.4 

285.1 ± 8.2 

282.1 

0.38 


TABLE XV. Experimental SO phase shift data used for the 
intermediate-energy polynomial fit, together with the fitted 
values and the residuals obtained from Eq. [S] The exper¬ 
imental data sets considered are the PY05 average [^, the 
re-analysis performed by Kaminski et al. [^, the CERN- 
Munich data of the Solution. B of Grayer et al. [sH] , and the 
results given by Hyams et al. [o^] . to which 5^°^ were added 
as systematic uncertainty. 
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Experiment 


(0“P ± 

^th (o) 

Ri 

Protopopescu 

1000 

0.37 ±0.08 

0.41 

-0.54 


1040 

0.38 ±0.04 

0.38 

-0.07 


1075 

0.45 ± 0.04 

0.48 

-0.67 


1105 

0.50 ±0.04 

0.55 

-1.36 


1135 

0.55 ± 0.04 

0.62 

-1.91 


1150 

.0.60 ±0.04 

0.66 

1.54 

Hyams 

969.2 

0.68 ±0.38 

1.00 

-0.83 


994.0 

0.08 ±0.08 

0.57 

-6.24 


1008.8 

0.08 ±0.14 

0.36 

0.15 


1028.6 

0.38 ±0.22 

0.36 

0.39 


1048.4 

0.61 ± 0.38 

0.40 

0.53 


1068.1 

0.84 ±0.18 

0.46 

2.14 


1087.9 

0.82 ± 0.21 

0.52 

1.43 


1107.7 

0.79 ± 0.22 

0.56 

1.04 


1127.5 

0.97 ±0.21 

0.61 

1.74 


1147.3 

0.88 ± 0.20 

0.66 

1.13 


1167.0 

0.87 ±0.22 

0.70 

0.77 


1186.8 

0.88 ± 0.20 

0.74 

0.72 


1206.6 

0.92 ± 0.21 

0.77 

0.71 


1226.4 

1.00 ± 0.21 

0.80 

0.94 


1246.2 

0.79 ± 0.15 

0.83 

-0.25 


1265.9 

0.84 ±0.14 

0.85 

-0.04 


1285.7 

0.95 ± 0.16 

0.86 

-0.52 


1305.5 

0.82 ± 0.16 

0.88 

-0.37 


1325.3 

0.87 ±0.16 

0.88 

-0.10 


1364.8 

0.82 ± 0.16 

0.89 

-0.47 


1384.6 

0.95 ± 0.17 

0.89 

0.34 


1404.4 

0.91 ± 0.25 

0.89 

0.08 


TABLE XVI. Experimental SO inelasticity data, together with 
the fitted values and the residuals obtained from Eq. [S] Only 
data sets from Hyams et al. and Protopopescu et al. 
which are consistent with the dip solution, are considered. 
With the exception of the data of Kaminski et al. as it is 
explained in the text above. 





